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利用 光 变 曲线 估计 猎 座 九 号 火箭 末 级 的 旋转 状态 研究 
AGER. WERE AXE UU, AgRuSRCC 
(1. 中 国 科 学 院 云南 天 文 台 ， 云 南 昆明 650011; 2. 中 国 科学 院 大 学 ， 北 京 100049; 3. 中 国 科学 院 空间 目 
标 与 碎片 观测 重点 实验 室 ， 江 苏 南京 210034) 

摘 要 : 光度 观测 是 地 基 观 测 空 间 目标 的 主要 观测 手段 之 一 ， 利 用 空间 目标 的 光度 信息 能 够 
估计 空间 目标 的 相关 特征 信息 . 为 了 更 好 地 了 解 空 间 目 标 旋 转 状态 , 选取 具有 代表 性 的 猫 座 
九 号 火箭 末 级 作为 研究 对 象 ， 由 其 光 变 信息 研究 旋转 状态 。 首 先 利 用 云南 天 文 台 1. 2 米 光学 
望远镜 获取 猫 座 九 号 火箭 末 级 的 光度 数据 ,再 对 目标 星 等 进行 斜 距 归 一 化 ,得 到 目标 光 变 信 
息 并 分 析 目 标 星 等 随时 间 变 化 的 光 变 曲线 ,估计 其 大 致 的 旋转 周期 ， 再 由 相位 离散 最 小 化 方 
法 计算 出 会 合 周期 。 根 据 太 阳 、 目 标 和 测 站 之 间 的 位 置 关系 、 惯 性 主轴 指向 、 旋 转轴 指向 、 
初始 相位 等 因素 ， 采 用 姿态 旋转 矩阵 来 计算 理论 星 等 ， 利 用 最 小 二 乘 原则 确定 惯性 主轴 方向 
及 初 相 角度 、 旋 转轴 指向 。 最 后 给 出 了 猎 座 九 号 火箭 末 级 的 旋转 周期 ,会 合 周期 以 及 旋转 轴 
指向 等 参数 ， 为 后 续 开展 其 他 空间 目标 光度 信息 研究 提供 参考 。 
关键 词 : 空间 碎片 的 旋转 状态 ; 光 变 曲线 ; 火箭 体 ; 相位 离散 最 小 化 ; 猫 座 九 号 ; 光度 模型 
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随 着 人 类 航天 活动 的 开展 ,地球 轨道 停留 了 越 来 越 多 的 空间 碎片 , 空间 碎片 的 存在 对 在 
轨 运 行 的 航天 器 以 及 载 人 航天 活动 带 来 很 大 的 威胁 。 其 中 火箭 末 级 也 属于 空间 碎片 , 是 指 将 
目标 航天 器 推 入 运行 轨道 的 最 后 一 级 火箭 , 在 完成 任务 后 会 绕 地 球 按 一 定 的 轨道 飞行 。 火 箭 
体 尺寸 、 重 量 大 且 运 行 姿态 不 受 控 制 , 会 严重 影响 附近 正常 工作 的 航天 器 ,一旦 相 接 会 产生 
更 多 的 空间 碎片 产生 连锁 效应 。 因 为 无 法 控制 火箭 体 的 姿态 和 轨道 , 为 了 了 解 这 些 大 型 空间 
碎片 的 运行 状态 , 需要 对 其 进行 观测 获取 来 其 相关 信息 。 地 基 光 度 观 测 是 最 常用 的 观测 手段 ， 
可 以 获取 目标 的 光度 信息 。 由 于 目标 的 光度 信息 与 目标 的 运动 、 姿 态 、 形 状 等 不 同 因素 县 奶 
相关 , 使 用 空间 态势 感知 的 不 同方 法 可 以 分 别 将 这 些 信息 分 离 出 来 ,从 而 得 到 目标 的 相关 信 
si 


IDo 


国内 外 分 析 空 间 目 标 光 度数 据 的 方法 主要 有 基于 相位 角 的 分 析 方法 、 基 于 极 值 点 的 分 析 
方法 ”、 基 于 光 变 曲线 周期 的 分 析 方 法 ”以 及 基于 模型 构建 的 分 析 方 法 等 。 前 三 种 方法 主 
要 根据 观测 得 到 的 光 变 曲线 进行 分 析 , 针对 不 同 种 类 的 空间 目标 分 析 方 法 也 有 所 不 同 。 但 是 
这 三 种 方法 的 适用 目标 具有 局 限 性 ,， 且 只 能 针对 某 一 特征 进行 分 析 。 基 于 模型 构建 的 分 析 方 
法 主要 包括 基于 非 线性 滤波 的 方法 ”、 基 于 封闭 面 元 模型 分 析 方法 和 基于 机 器 学 习 方法 等 ， 
其 中 基于 非 线 性 滤波 的 方法 是 分 析 空 间 目 标 特性 的 主要 研究 方法 , 通过 构建 目标 运动 和 观测 
模型 来 预测 其 运动 和 姿态 , 需要 精确 的 先 验 参数 才能 得 到 收敛 的 结果 , 不 适用 于 先 验 参数 较 
少 的 目标 。 本 文 主要 根据 观测 到 的 火箭 体 光 度数 据 , 对 处 理 后 的 光 变 曲线 进行 周期 分 析 来 确 
定 火 箭 体 的 旋转 周期 , 通过 目标 相位 角 和 形状 参数 建立 光度 模型 , 通过 比 对 观测 光 变 曲线 和 
蛙 论 计算 得 到 的 光 变 曲线 估计 旋转 轴 的 指向 ,根据 得 到 的 收敛 区 间 证 明了 构建 光度 模型 方法 
的 可 行 性 ， 为 后 续 研 究 其 他 类 型 的 空间 目标 姿态 研究 提供 参考 。 

* 基金 项 目 ， 国家 自然 科学 天 文联 合 基金 〈U1731112，U1831106) ; 国家 自然 科学 基金 (11655001，11703011〉; 中 国 科 学 院 
青年 创新 促进 会 项 目 资助 (Y201817〉; 云南 省 自然 科学 基金 重点 项 目 (2019FA002) . 
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普 鹰 九 号 火箭 末 级 〈 编 号 ; 429680 作为 研究 目标 。 猪 鹰 九 号 火箭 是 SpaceX 


公司 研制 的 重型 运载 火箭 , 是 现役 运载 能 力 最 大 的 火箭 之 一 ， 且 其 末 级 火箭 的 太 寸 较 其 他 火 
第 体 更 大 , 更 易 进 行 预报 和 进行 光度 观测 ,同时 大 尺寸 也 使 得 其 对 于 正常 工作 的 航天 器 影响 


更 大 。 了 解 其 姿态 可 以 更 好 地 预测 它 的 运动 从 而 使 其 他 正常 工作 的 航天 器 更 好 地 规避 碰撞 风 
险 。 截止 到 20204£ 5 H 30 日 已 经 发 射 了 83 次。 根据 SpaceX 公司 发 布 的 数据 ，42968 =H 


标 是 猎 鹰 九 号 中 v1.2 型 火箭 的 火箭 末 级 , 如 图 1 所 示 。 其 去 除 下 底面 引擎 后 长 约 为 9. 44m， 
底面 直径 为 3. 66m。 其 轨道 为 大 椭圆 轨道 ， 近 地 点 230km， 远 地 点 33819km。 有 具体 信息 如 表 1 
所 示 : 
表 1 目标 42968 相关 信息 
Table 1 informations of target 42968 
目标 名 称 国际 编号 类 型 发 射 日 其 轨道 周 其 轨道 倾角 近地点 远地点 
FALCON 9 R/B 2017-063B 火箭 体 2017-10-11 593. 83min 27.80 230km 33819km 
FALCON 9 V1.1 FALCON 9 V1.2 
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图 1 猎 鹰 九 号 火 科 尺寸 信息 


Fig 1 size information of Falcon 9 Rocket 


2. 光度 数据 的 获取 流程 


本 文采 用 地 基 光 度 测 量 的 方法 获取 目标 的 光度 信息 ,利用 云南 天 文 台 1.2m 望远镜 对 目标 
进行 观测 。 将 采集 到 的 图 像 进行 本 底 和 平 场 改 正 、 流 量 定 标 、 斜 距 归 一 化 等 处 理 后 得 到 目标 


的 视 星 等 及 星 等 随时 间 变 化 的 光 变 曲线 。 
2.1 图 像 采集 
观测 前 从 space-track 网 站 上 获取 卫星 的 轨道 数据 ,一 般 为 双 行 根 数 格式 (TLE)。TLE 的 


预报 精度 不 高 ， 应 用 到 激光 测 距 中 需要 修正 ， 张 训 方 等 人 研究 了 实时 修正 预报 的 方法 05] ， 
而 本 文 观测 使 用 的 望远镜 视 场 较 大 ， 可 以 不 进行 轨道 修正 。 使 用 云南 天 文 台 1.2m 望远镜 及 
30cm 导 星 镜 、 配 套 控 制 系统 、CCD 相机 对 目标 进行 观测 并 采集 平 场 、 定 标 星光 度 信 息 及 
目标 的 光度 信息 ， 使 用 SourceExtractorr4 软 件 提取 定 标 星 和 待 测 目标 的 流量 值 信息 。 
观测 时 间 : 2019 年 12 月 13 A , 2019 Ẹ 12 月 15 日 ,每 次 观测 时 间 约 为 30 分 钟 ， 具 
体 信 息 如 表 2 所 示 。 图 2 是 使 用 1.2m 望远镜 拍摄 的 目标 光度 图 像 ,处 于 图 像 中 间 的 亮点 即 
为 观测 目标 。 图 3 表示 13 日 和 15 日 观测 时 段 内 目标 相对 于 望远镜 方位 角 和 高 度 角 随时 间 
变化 的 曲线 ， 以 极 坐 标的 形式 表示 出 来 。13 日 用 实 线 表 示 ，15 日 用 虚线 表示 ,每 段 曲线 的 


初始 方位 和 俯仰 角 通 过 曲线 中 标注 的 点 表示 出 来 ,可 以 看 出 在 观测 时 段 内 目标 的 方位 和 高 度 


变化 并 不 大 。 
Tab.2 the information of observation 
观测 日 期 观测 距离 观测 时 间 时 长 曝光 时 间 
2019 年 12 A 13 A 20458. 54—24682. 70 12:21:59. 86——12:52:23. 33 31min 0. 5s 


2019 ££ 12 H 15 H 17166. 77——21655. 37 14:08:25. 15-—-14:37:33. 19 29min 0. 5s 
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图 2 目标 光度 图 像 


Fig2 object’s luminosity figure 


The elevation and azimuth of 42968 in observe time 
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图 3 观测 时 段 内 目标 方位 角 与 俯仰 


Fig3 space object’s azimuth and elevation in observe periods 


2.2 平 场 、 本 底 处 理 及 流量 定 标 


在 采集 到 本 底 与 平 场 图 像 后 进行 归 一 化 处 理 。 计 算 经 过 平 场 处 理 后 的 图 像 厌 度 值 公式 如 


其 中 6o 是 原始 图 像 的 厌 度 值 ，Gpics 是 本 底 厌 度 值 ，Grrut 是 标准 平 场 厌 度 值 
根据 现 有 的 星 表 选择 定 标 星 并 进行 观测 ， 由 于 望远镜 视 场 较 小 ， 不 能 用 场 星 定 标 ， 因 此 
本 文选 择 的 是 landolt 星 表 里 的 恒星 G162 , 采集 到 其 光度 图 像 后 经 过 本 底 、 平 场 处 理 得 到 


处 理 后 的 图 像 ,这 一 过 程 是 消除 望远镜 光学 系统 对 于 目标 光度 的 影响 。 将 处 理 后 的 图 像 与 星 


表 中 的 位 置 进行 比 对 来 确定 恒星 ， 得 到 恒星 光度 流量 值 ,再 除 以 其 曝光 时 间 ,得 到 定 标 星 单 
位 时 间 内 的 流量 值 六 。，, 其 视 星 等 值 已 知 为 Muw=13.012 ,这 样 即 可 以 获得 仪器 的 星 等 替 点 : 
Mo = 2.5: 10810 (fcar) + Mcai (2) 
从 而 可 以 计算 得 到 目标 的 视 星 等 : 
Mopj = —2.5: logio(fop;) + Mo (3) 


其 中 记 w) 为 待 测 目标 的 流量 值 。 


2. 3 异常 点 及 斜 距 归 一 化 处 理 


在 使 用 SourceExtractor 软件 提取 流量 过 程 中 会 识别 到 非 待 测 目标 的 流量 信息 ， 从 而 得 
到 错误 的 星 等 值 ,在 观测 过 程 中 望远镜 对 目标 进行 跟踪 ,因此 目标 在 图 像 中 的 位 置 相 对 固定 ， 
因 预 报 误差 而 缓慢 变化 ,因此 可 在 目标 轨道 预报 精确 或 目标 运动 速度 不 快 的 情况 下 可 以 根据 
目标 在 fits 图 像 中 的 坐标 将 错误 识别 的 目标 剔除 如 图 4 所 示 ， 从 而 计算 得 到 正确 的 星 等 。 
4 表示 了 13 日 及 15 日 的 异常 点 在 CCD 坐标 的 相关 信息 ， 可 以 看 到 一 部 分 点 离散 在 目标 坐 
标 区 域外 ,根据 这 一 信息 将 坐标 错误 的 点 剔除 。 这 一 方法 不 能 将 被 恒星 污染 、 云 层 遮 挡 这 样 
的 情况 剔除 掉 ， 而 这 两 种 情况 下 目标 的 星 等 值 会 发 生 突变 ， 所 以 可 以 将 距离 平均 光 变 曲线 
0.5 星 等 的 数据 点 fits 图 像 与 在 平均 光 变 曲线 上 数据 点 的 fits 图 像 进行 比 对 ， 确 认得 到 的 是 
否 是 待 测 目标 的 光度 信息 图 4 中 15 日 目标 坐标 区 域 中 的 萎 形 点 就 是 被 恒星 污染 的 异常 点 。 

将 异常 点 剔除 后 得 到 了 关于 目标 的 光度 数据 。 为 了 去 除 待 测 目标 斜 距 对 于 目标 光度 的 影 
响 ， 需 要 对 得 到 的 星 等 进行 斜 距 归 一 化 处 理 ， 通常 将 卫星 斜 距 归 一 化 到 1000km 处 。 根 据 


公式 : 


Rob; 
Muni = Mopj — 5-logio( 9) (4) 


Runi 
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得 到 归 一 化 后 的 星 等 Mn;。 其 中 Robj 为 卫星 距 测 站 的 斜 距 ， Runi = 1000km。 

5 为 剔除 异常 点 、 斜 距 归 一 化 后 的 光 变 曲线 ,其 中 横 坐 标 表 示 时 间 ,以 s 为 单位 ; 纵 
坐标 表示 星 等 。 由 图 可 见 光度 随时 间 的 变化 呈现 出 周期 性 特征 ,而且 在 观测 时 段 内 星 等 的 变 
化 周期 基本 保持 不 变 ， 大 致 为 125s。 整 体 光 变 曲线 曲线 也 有 一 个 缓慢 的 变化 趋势 ， 这 两 天 
的 星 等 整体 有 缓慢 增 大 的 趋势 ， 曲 线 振幅 在 逐渐 变 大 。 下 一 步 会 针对 这 两 条 光 变 曲线 进行 分 


析 以 确定 目标 的 旋转 速率 及 旋转 轴 的 指向 。 


Coordinate of Outlier and Normal Point 
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图 4 目标 异常 点 坐标 分 布 


Fig4 Outlier of space object's coordinate distribution 


Lightcurve of Satellite 42968 in 20191213 


-0.2 ii Lightcurve of Satellite 42968 in 20191215 
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图 5 观测 时 段 内 目标 42968 的 光 变 曲线 


Fig.5 The lightcurve of 42968 in observation periods 


3、 火 箭 体 旋转 状态 分 析 
火箭 体 属 于 空间 碎片 且 其 运行 轨道 和 运行 姿态 不 受 地 面 控制 。 描 述 火箭 体 旋转 状态 主要 包 
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括 旋 转 周期 、 旋 转轴 指向 和 旋转 初 相位 这 几 个 参数 。 本 节 主 要 根据 第 2 节 得 到 的 目标 光 变 


曲线 来 确定 目标 的 旋转 周期 和 旋转 轴 指 向 ， 并 分 析 得 到 的 结果 。 


3.1 会 合 周 期 分 析 

会 合 周期 是 由 一 种 相对 周期 , 取决 于 卫星 、 测 站 及 太阳 的 位 置 关 系 。 根据 网 球拍 定理 
可 知 , 自由 刚体 会 逐渐 绕 着 目标 最 大 转动 惯量 的 轴 旋 转 。 42968 号 火箭 末 级 入 轨 时 间 较 长 (从 
2017 年 10 月 发 射 至 观测 时 刻 约 两 年 )， 其 几何 形状 近似 为 圆柱 体 ， 所 以 可 以 合理 假设 目标 
绕 着 垂直 于 目标 侧面 的 轴 旋 转 ， 对 于 测 站 来 说 观测 到 的 星 等 变化 周期 上 只是 旋转 周期 的 一 半 ， 
故 其 会 合 周期 大 约 是 250s。 为 了 更 加 精确 地 确定 会 合 周期 ， 本 文采 用 了 相位 离散 最 小 化 方 
法 (Phase Dispersion Minimization, 简称 PDM) “将 数据 点 分 段 进行 分 析 得 到 最 优 解 。PDM 
的 原理 是 将 观测 数据 按 试 验 周期 计算 , 根据 这 些 试验 周期 的 计算 相位 矢量 , 将 这 些 相 位 矢量 
按 (0,1) 区 间 划 分 为 不 同 的 子 区 间 ， 计 算 子 区 间 与 整个 区 间 的 数据 方差 ， 使 得 二 者 之 比 最 小 
的 试验 周期 即 是 最 精确 的 会 合 周期 。 

空间 目标 的 离散 观测 值 可 以 分 为 两 个 矢量 : 星 等 矢量 风 及 观测 时 间 矢 量 #， 其 中 第 i 个 
观测 点 可 以 表示 为 :Cm ti), 假 设 在 一 个 观测 时 段 内 共有 NN 个 观测 离散 点 , 即 i = 1,2,...,N。 则 
星 等 矢量 元 的 方差 ?为 : 


2 ..Xmi-m) 
EDS (5) 


Jem = O9 mft fl 
在 星 等 矢量 元 中 选取 个 不 同 的 数据 子 集 , FL UEREAM BOR Ten p Grm 4 Hi EUER 
长 的 样本 方差 为 s? (j=1, 2,…, MD 。 则 所 有 数据 集 的 总 方差 s? 为 : 


X(n;-1)s? 
2. p og 
dao 6) 


in 


> 


假设 BP 为 试验 周期 矢量 ，B 中 的 元 素 为 pe， 根 据 试验 周期 计算 时 间 矢 量 ? 的 相位 矢量 gp。 
对 于 相位 矢量 2 中 的 每 个 元 素 内 来 说 ， 


= A= = (7) 
o PERRERA BAR qi A H E i E v E HG £s 


假设 从 元 中 选取 的 M 个 子 数据 集 有 相似 的 相位 矢量 Bp。 将 全 相位 间隔 (0, 1) 分 为 固定 大 
小 的 等 长 区 间 。 定 义 统计 量 : e= 三 ， 其 中 ga? 和 s? 分 别 由 公式 6) 和 公式 (7) 得 到 。 如 果 
px 个 是 正确 的 周期 ，8 = 1; 如 果 px 是 正确 的 周期 ，8 在 试验 周期 内 会 达到 最 小 值 ， 晶 趋 近 
FO, 

由 火箭 体 光度 数据 ，t 为 观测 时 段 的 时 间 点 ， 元 为 观测 时 段 的 目标 星 等 值 ， 会 合 周期 大 
约 为 250s， 所 以 试验 周期 矢量 B 取 值 为 200 一 300， 步 长 取 0. 1。 相 位 间隔 步 长 取 0.1， 即 将 
(0，1) 区 间 划 分 为 10 个 子 区 间 。 使 用 上 述 参数 及 相位 离散 最 小 化 方法 计算 得 到 更 为 准确 
的 火箭 体会 合 周 期 。 

图 6 表示 了 两 个 观测 时 段 内 试验 周期 及 通过 相位 离散 最 小 化 方法 计算 得 到 的 试验 周期 
对 应 的 9 值 。 由 图 6 的 两 条 曲线 可 以 发 现在 试验 周期 内 都 存在 一 个 极 小 值 点 ，13 日 在 6 最 小 
处 时 对 应 的 周期 为 249. 23s, 由 图 6 中 的 下 图 可 以 得 到 15 日 在 9 最 小 处 对 应 的 周期 为 248. 29s。 
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图 6 火箭 体会 合 周 期 分 析 


Fig.6 Analysis of Rocket Body's Synodic Period 


3.2 旋转 轴 指 向 分 析 


空间 目标 的 旋转 状态 包含 旋转 周期 及 旋转 轴 指 向 ， 其 会 合 周期 已 经 通过 PDM 方法 得 到 ， 
对 于 此 观测 目标 , 其 星 等 变化 周期 比较 稳定 且 相 对 于 测 站 及 太阳 的 相位 角 变 化 不 大 , 所 以 该 
目标 的 会 合 周期 可 近似 看 作 其 旋转 周期 。 

在 赤道 惯性 坐标 系 下 旋转 轴 以 矢量 BP 表示， 转化 为 赤 经 g、 赤 纬 6 表 示 : 


Pxi cosócosa 
P-|P*yi| 2|cosósina 
P*z sinó 


在 分 析 旋 转轴 指向 之 前 需要 建立 目标 光度 模型 。 火 箭 体 本 身 并 不 发 光 而 是 反射 太阳 光 ， 
测 站 观测 到 的 目标 图 像 其 实 是 火箭 体 的 反射 光 ， 与 目标 表面 材料 、 目 标 形 状 等 因素 相关 、 太 
阳 一 目标 一 测 站 之 间 的 几何 关系 相关 ， 三 者 的 几何 关系 如 图 7 所 示 ， 其 中 wps 和 Usun 分 别 表 
示 目 标 到 地 球 和 太阳 的 方向 单位 矢量 ， 友 表示 目标 表面 面 元 的 法 线 方向 单位 矢量 ，Uoms 和 
Usum 之 间 的 夹 角 即 相位 角 ipy, 玉 是 相位 角 平 分 线 方向 单位 矢量 。 由 目标 的 双 行 根 数 得 到 其 位 
置信 息 ，13 日 的 双 行 根 数 历 元 为 2019-12-12T18:50:17Z2，15 日 的 双 行 根 数 历 元 为 
2019-12-14T10:29:092， 太 阳 的 位 置信 息 由 JPL 的 de430 星 历 给 出 ， 测 站 位 置 : (25.0299 
N, 102. 7974 E) ， 海 拔高 度 为 1991. 83m。 将 这 三 者 的 位 置 转换 到 惯性 系 下 表示 为 位 置 矢 量 ， 
根据 三 者 的 位 置 矢量 即 可 求 出 Wops，Usun， 则 相位 角 ipy: 


" u, 1 
ip, = arccos( —5 —u- (8) 
[uopslIusunl 


观测 时 段 内 相位 角 的 变化 如 图 8 所 示 . 


二 
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Satellite Surface 


图 7 太阳 -目标 - 测 站 之 间 的 儿 何 关 系 


Fig.7 The geometry relationship of sun-space object-observation station 


Phase Angle of Satellite 42968 
Ue 20191213 
—— 20191215 
24 


N 
N 


N 
o 


Phase Angle of Satellite(deg) 


m 
[] 


45000 46000 47000 48000 49000 50000 51000 52000 53000 
Second of Day(s} 


图 8 观测 时 段 内 目标 相位 角 变 化 曲线 


Fig.8 variation trend of object's phase angle in observation period 


火箭 体 的 形状 可 看 作 圆 柱 体 , 根据 此 假设 可 以 计算 火箭 体 的 等 效 面 积 。 假 设 太 阳光 在 火 
第 体 表 面 为 漫 反 射 ， 圆柱 体 侧面 的 法 线 方向 定义 为 : 圆柱 的 上 下 底 中 心 连 线 ， 向 上 底 方向 为 
1E; 圆柱 体 侧面 面 元 的 法 线 方向 为 垂直 于 面 元 方向 向 外 , 对 圆柱 体 侧面 面 元 积分 就 可 以 得 到 
侧面 的 等 效 面积 。 大 阳 入 射 角 为 isun， 卫 星 一 测 站 方向 矢量 与 法 线 方 向 矢量 夹 角 为 iops。 则 
圆柱 体 侧 面 等 效 效 面积 关 


Acyea = | sinis) coso ' sin(igp,)cos(q — 0) ‘hsdg = EE [sing + (1 — 
s 
0): cos6] (9) 
其 中 0 满足 : 
eosdem Cos(ipa)—cos(isun)*cos(iobs) (1 0) 
sin(isun)*sin(iobs) 


HPA A REW TERR, oA EAER MAEA, hA REAA, dA E 
体 底面 直径 。 
圆柱 体 上 下 两 底 可 近似 看 作 平面 ， 平 面 的 等 效 面 积 为 : 
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Aplea Api “cos(isun)cos(iops) (11) 


不 考虑 下 底 的 火箭 体 喷 嘴 ,， 火 箭 体 上 下 底 的 法 线 方向 分 别 沿 着 上 下 底 平 面 从 里 向 外 ,分别 使 
用 公式 (11) 计算 火箭 体 上 下 底 的 等 效 面 积 。 根 据 火 箭 体 侧面 及 上 下 底 的 等 效 面积 就 可 以 确定 
火箭 体 的 等 效 面积 从 而 计算 目标 的 模型 光度 。 

给 定 火箭 体 姿 态 初始 参数 : Ao，Ai : 目标 惯性 主轴 与 旋转 轴 夹 角 ; Ao. 09: 目标 旋转 轴 
指向 ，Xo 即 惯性 系 的 赤 经 ，00 为 惯性 系 的 余 纬 (旋转 轴 与 惯性 系 z ef o: 旋转 初 
始 相位 ; To: 目标 旋转 周期 。 随 着 时 间 的 变化 相位 8 也 在 变化 ，% = po+at. 至 ，dt 表示 时 间 的 
变化 。 将 Ao, AM, Ao, 0o, 通过 姿态 旋转 矩阵 从 惯性 坐标 系 转换 到 星体 坐标 系 下 表示 ， 并 将 这 
些 旋转 矩阵 相 乘 得 到 新 的 和 矩阵 R，R 和 表面 法 向 量 的 点 积 即 为 不 同时 刻下 法 向 量 在 星体 坐标 
系 下 的 指向 。 再 根据 太阳 单位 入 射 辐射 流量 密度 照射 在 有 效 面积 上 反射 的 流量 及 太阳 的 视 星 
等 计算 出 模型 星 等 值 。 根据 前 文 所 述 ,， 旋转 周期 近似 采用 会 合 周 期 ,目标 绕 最 大 转动 惯量 轴 
( 短 轴 ) 转动 ， 即 惯性 主轴 与 旋转 轴 夹 角 Ao= 0, A1= 90, 旋转 周期 7。 = 249.23s 只 需要 确 钢 
述 旋转 轴 指 向 的 两 个 角 4o，6o 及 旋转 初 相 bo， 通 过 参数 搜索 ， 以 最 小 二 乘 为 判 据 确 定 旋转 
参数 。 

以 13 日 的 数据 为 例 ,已 知 Ao= 0 ,Ai= 90 ,To = 249.23s。 遍 历 三 个 参数 :10 € [0,360], 
步 长 为 2; 6o € [0,180]， 步 长 为 1; ġo € [0,360]， 步 长 为 10。 每 次 循环 计算 观测 时 段 内 观 
测 星 等 与 模型 星 等 之 差 的 方差 。 裔 历 完 所 有 参数 后 得 到 总 体 方差 分 布 ， 以 Bo 划分 共有 36 个 
0-C 方差 分 布 图 , 在 每 个 方差 分 布 图 中 横 坐 标 表 示 旋 转轴 指向 4o; 纵 坐 标 表 示 旋 转轴 指向 06， 
发 现 每 张 0-C 方差 分 布 图 都 有 两 个 方差 极 小 值 区 域 , 分 别 位 于 每 张 子 图 中 的 左下 方 及 右上 方 ， 
且 关 于 中 心 点 对 称 。 这 是 由 于 建立 光度 模型 时 假设 火箭 体 上 底 及 下 底面 积 相等 导致 的 对 称 性 ， 
根据 这 一 推论 可 知 在 两 个 方差 极 小 值 区 域 分 别 存在 使 得 方差 最 小 的 参数 组 合 。 在 每 个 方差 极 
小 值 区 域内 都 有 两 个 最 小 方差 收敛 区 域 ，0-C 方差 最 小 值 点 在 其 中 一 个 方差 收敛 区 间 内 。 

计算 每 张 等 高 线 图 内 最 小 方差 及 所 对 应 的 (Xo0,00,80) 组 合 ， 得 到 了 36 个 最 小 方差 点 ， 
再 在 这 36 种 组 合 中 选择 最 小 方差 对 应 的 组 合 : 9,09, bo = (212,116,120)。 图 9 中 左 图 给 出 
了 最 小 方差 组 合 所 对 应 的 方差 分 布 图 , 等 高 线 颜 色 越 接近 蓝 色 表示 方差 越 小 , 越 接近 红色 表 
示 方 差 越 大 。 根 据 0-C 方差 的 对 称 性 可 知 在 方差 分 布 图 的 左下 方 也 应 存在 使 得 0-C 方差 最 小 
的 参数 组 合 ， 计 算 关 于 左下 方差 极 小 值 区 域 0-C 方差 最 小 值 的 参数 组 合 ， 得 到 了 图 9 AR, 
方差 极 小 值 点 对 应 的 参数 为 ，40,00, Po = (32,64,290)， 处 于 0-C 方差 收敛 区 间 内 。 

根据 以 上 得 到 的 两 个 参数 组 合 , 将 其 分 别 代 入 到 光度 模型 中 , 并 比 对 模型 光 变 曲线 和 观 
测 光 变 曲线 ， 比 对 结果 如 图 10 所 示 。 图 10 中 的 左右 图 分 别 对 应 图 9 中 左右 图 得 到 的 参数 ， 
图 中 标识 出 了 每 个 参数 组 合 对 应 的 0-C 方差 值 , 可 以 发 现 左 图 0-C 方差 值 更 小 。 不 过 右 图 的 
方差 与 左 图 方差 相差 并 不 大 并 且 目 标 星 等 的 变化 趋势 也 很 好 的 拟 合 了 出 来 ,以 左 图 为 例 可 以 
看 到 光度 模型 将 火箭 体 光 度 整体 的 变化 趋势 表示 了 出 来 ,包括 振幅 逐渐 变 大 的 趋势 及 每 个 周 
期 星 等 极 大 值 点 的 变化 .不 过 对 于 每 个 周期 星 等 极 小 值 点 处 的 光度 突变 光度 模型 并 没有 很 好 
的 拟 合 出 来 , 这 一 部 分 的 突变 可 能 是 由 于 火箭 体 侧面 镜 反射 导致 的 ， 而 光度 模型 上 只 考虑 了 目 
标 漫 反射 的 情况 。 

根据 以 上 分 析 可 以 确定 卫星 编号 为 42968 的 目标 旋转 轴 指 向 可 能 包含 两 种 组 合 , 分 别 是 
Ao, 09 = (212,116) 及 1o, 09 = (32,64)， 对 应 的 旋转 初 相 位 分 别 是 120，290， 二 者 之 差 为 170 
度 ， 而 且 求 得 的 两 个 旋转 轴 指 向 赤 纬 416 相 差 180 度 ， 两 个 余 纬 90 与 90” 之 差 的 绝对 值 都 为 
26 度 ， 正 好 对 应 了 光度 建 模 时 火箭 末 级 上 下 底面 积 相 等 而 导致 求 得 的 两 个 旋转 轴 指 向 关于 
惯性 系 原点 对 称 的 情况 。 
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图 10 13 日 模型 及 观测 光 变 曲线 


Fig.10 lightcurve of observation and calculate in 2019.12.13 


本 文通 过 分 析 猪 鹰 九 号 火箭 体 (NORAD 编号 : 42968) 的 光 变 曲线 确定 其 旋转 状态 。 根 据 
2019 年 12 月 13 日 和 15 日 实测 的 数据 得 到 的 目标 光 变 曲线 存在 明显 的 周期 性 。 这 说 明 目 标 
的 运动 姿态 较为 稳定 。 由 于 在 两 天 的 观测 时 段 内 目标 相位 角 的 不 同 造成 了 归 一 化 处 理 后 的 目 
标 绝对 星 等 的 细微 差别 。 由 光 变 曲线 的 周期 并 通过 PDM 确定 了 目标 的 会 合 周期 : 13 日 该 目 
标的 会 合 周期 为 249. 23s，15 日 该 目标 的 会 合 周期 为 248. 29s， 在 两 天 的 时 间 里 目标 的 周期 
并 没有 发 生 的 明显 变化 , 通过 建立 光度 模型 的 方法 计算 目标 的 模型 星 等 , 搜寻 旋转 轴 指 向 和 
旋转 初始 相位 ,以 模型 星 等 与 观测 星 等 之 差 的 方差 最 小 为 判 据 确定 目标 的 旋转 轴 指 向 。 本文 
计算 了 13 日 的 / C 方差 ， 由 于 模型 假设 火箭 体 的 上 底 及 下 底面 积 相等 所 以 得 到 了 两 个 可 能 
的 旋转 轴 指 向 ， 分 别 为 1o, 09 = (212,116) 及 4o, 09 = (32,.64)， 若 要 更 加 精确 地 确定 旋转 轴 指 
向 还 需 更 多 的 观测 数据 及 完善 光度 模型 。 
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Abstract: Lightcurve measurement is one of the main methods for ground-based space-target 
Observation, as the photometric information of the targets can be used to estimate their relevant 
characteristics. Aiming at better understanding of the rotational state of the space target, a 
representative upper stage of the Falcon-9 rocket is selected as the object to analyze its rotational 
state through light curve study. Firstly, the photometric data of the upper stage of the Falcon-9 
rocket is obtained with the 1.2-m optical telescope of the Yunnan Observatories. Then the target's 
magnitude is processed with slant-range normalization. The photometric changing information is 
obtained, and the time-dependent light curve of the target magnitude is analyzed. Accordingly, the 
approximate rotation period of the target is estimated, and then the accurate synodic period is 
calculated by the phase discrete minimization method. Based on the position relationship between 
the Sun, the target and the observatory, the main inertial axis pointing, the rotation axis pointing, 
the initial phase and other factors, the theoretical magnitude is calculated using the attitude 
rotation matrix. The main inertial axis direction, the initial phase angle and the rotation axis 
pointing direction are determined as well. Finally, the parameters such as the rotational period, the 
synodic period and the pointing direction of the rotational axis of the upper stage of the Falcon 9 
rocket are given, which can serve as a reference for following studies on the photometric 


information of other space targets. 
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